Exploring the fragile antiferromagnetic superconducting phase in CeCoIn 5 
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CeCoIns is a heavy fermion Type-II superconductor which exhibits clear indications of Pauli- 
limited superconductivity. A variety of measurements give evidence for a transition at high magnetic 
fields inside the superconducting state, when the field is applied either parallel to or perpendicular 
to c axis. When the field is perpendicular to the c axis, antiferromagnetic order is observed on the 
high- field side of the transition, with a magnetic wavevector of (q q 0.5), where q = 0.44 reciprocal 
lattice units. We show that this order remains as the magnetic field is rotated out of the basal 
plane, but the associated moment eventually disappears above 17°, indicating that the anomalies 
seen with the field parallel to the c axis are not related to this magnetic order. We discuss the 
implications of this finding. 

PACS numbers: 74.70.Tx,75.30.Kz 



CeCoIns is one of the more prominent members of the 
Ce family of heavy fermion compounds. It is an ambi- 
ent pressure superconductor with a zero field T c of 2.3 K, 
the highest of any heavy fermion material [l[ . It is an ex- 
tremely clean high-K superconductor, with an electronic 
mean free path ~ 1 ^m @, It is a strongly Type- 
II superconductor, and there is clear evidence that the 
superconductivity is Pauli-limited in strong fields 

The material is anisotropic; it has the tetragonal 
HoCoGa 5 structure, with a = 4.602 A and c = 7.545 A at 
2 K (space group PA/mmm). If the field is applied in 
the basal plane, the zero-temperature upper critical field 
HoH C 2 — 11.6 T, but when the field is parallel to the 
c axis, ilqH C 2 = 4.95 T. However, in both cases, it en- 
ters the Pauli-limited region close to H C 2 0. A vari- 
ety of measurements at high fields have found evidence 
for a transition inside the superconducting mixed state 
for both field orientations. These measurements include 
specific heat @,J3, ultrasound studies of the transverse 
shear velocity [J, thermal conductivity NMR 
and magnetostriction [Io| . This transition, existing in 
a high-field, low-temperature region bordering H C 2 has 
been considered as a candidate for a Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state [HI, [Hj]. In this state, the 




Zeeman splitting of the energies of electronic quasiparti- 
cles is expected to lead to an additional spatial variation 
of the superconducting order parameter, for instance [l2[ 
with nodal planes perpendicular to the flux lines. 

Kenzelmann et al. (i~3l . fbl ] have shown that when a 
magnetic field is applied in the basal plane there is in- 
commensurate antiferromagnetic order inside the high- 
field phase. Magnetic Bragg reflections are seen with a 
wavevector ±(q q 0.5) relative to a nuclear peak where 
q = 0.44 reciprocal lattice units (rlu) when the field is 



Figure 1: A sketch of the geometrical constraints imposed 
by the cryomagnet. ki is the incident neutron wavevector 
and k/ is the outgoing neutron wavevector. For this sample 
geometry, absorption is minimised if the outgoing (rather than 
the incoming) beam is nearly parallel to the sample plate. The 
circular arrow indicates the direction of tilt required to rotate 
the magnetic reflections into the scattering plane. 



applied parallel to [110]. By symmetry, the wavevector 
is (q -q 0.5) when the field is parallel to [1 1 0]. The 
magnetism is carried by the Ce ions, and was found to 
have the same value of 0.15 ± 0.05 /xs/Ce ion and the 
same ordering vector for two different field directions in 
the basal plane. This antiferromagnetic order disappears 
sharply with increasing field as the superconductivity is 
destroyed at H c2 . If the magnetic field is applied along [1 
0], the two wavevectors are degenerate and the sample 
is expected to have domains of both [lij]. An alterna- 
tive scenario of double-q magnetic ordering is ruled out 
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Figure 2: (color online) The scattering plane with the cry- 
omagnet level. The thick dashed lines show the two field 
orientations examined in this paper. The marked areas are 
the reciprocal space accessible with a wavevector of 1.4 A" 1 
(the wavevector used to gather the data presented here). The 
open circles represent nuclear reflections and the closed points 
represent magnetic reflections, although the magnet had to be 
tilted to bring them into the scattering plane (see main text). 
For measurements of the nuclear reflections, the direction of 
the field is irrelevant, and so the sample could be rotated 
freely to bring them into the window of accessible reciprocal 
space. 



by the observation of a characteristic single-q NMR line- 
shape [l5j ]. 

We have carried out neutron diffraction measurements 
on a sample of CeCoIns, to investigate changes in the 
antifcrromagnetic order when the applied field is rotated 
out of the basal plane. Our scattering geometry also 
allows the coherence of the order along the field direc- 
tion to be investigated. The measurements reported here 
cast further light on the relationship between the mag- 
netic and superconducting order parameters in CeCoIiis, 
which remains a matter of debate. 

To achieve this geometry, we required a horizontal 
magnetic field, so that it could be in the scattering plane, 
and of sufficiently large magnitude to probe the region of 
interest. For this reason, our experiment could only be 
carried out at the RITA-II instrument at the SINQ Facil- 
ity (Paul Schcrrcr Institut), using the 11 T split-pair hor- 
izontal field cryomagnet available there, equipped with a 
dilution refrigerator (DR) insert. This magnet has two 
neutron access ports with an acceptance angle of ± 15°, 
parallel to the field, and two with an acceptance angle of 
± 7.5°, perpendicular to the field (Figure [T]). The mag- 
net could be rotated about the vertical axis and tilted a 
few degrees about a horizontal axis. The DR could be 
rotated independently about the vertical magnet axis to 
alter the orientation of the sample relative to the field 



(lj). However, the limited neutron access (inherent in 
any high field magnet) imposed the constraint of 29 = 
90° ± 22.5° scattering, severely limiting the accessible 
region of reciprocal space (see Figure [21). 

CeCoIn5 crystals grow preferentially as thin plates 
with the c axis normal to the plate. Our 146 mg sin- 
gle crystal was glued to an aluminum plate mounted on 
the DR to give a [1 1 0]-[0 1] horizontal scattering plane 
when the cryomagnet was level; the DR was rotated to 
bring the field direction close to [1 1 0] . Figure [2] shows 
the accessible regions of reciprocal space in this scattering 
plane for two orientations of the sample relative to the 
field. A typical magnetic reflection close to this plane is 
at (1-q, -q, 1.5) = (0.56 -0.44 1.5). To bring the magnetic 
reflections into the scattering plane, the cryomagnet was 
tilted by a few degrees about the horizontal axis perpen- 
dicular to the field. 

RITA-II is a cold three axis spectrometer, with a set of 
seven analyzer crystals ('blades') receiving the scattered 
beam These blades reflect into a 2-D detector, and 
can be set up in a variety of arrangements. We operated 
the instrument in its elastic scattering configuration [l6| , 
using the blades to probe multiple 29 values; we assume 
that all scattering observed is elastic in origin. The spec- 
trometer and sample were always aligned such that the 
Bragg reflection being studied was incident on the cen- 
tral blade. A beryllium filter and radial collimator were 
installed on the output beam. Blade efficiency was cor- 
rected for by using incoherent scattering gathered under 
conditions of field or blade angle for which no magnetic 
signal was observed. 

Figure [3] shows the magnetic reflection (0.56 -0.44 1.5) 
measured at T = 50 mK as a function of the sample 
rotation angle lo with a field of 10.8 T applied out of 
the basal plane at 11.9° to [110], known with an accu- 
racy of ± 0.1°. The background here is a composite of 
data obtained at neighbouring 29 values at 10.8 T and at 
10.9 T. The peak has been fitted as a Gaussian, giving a 
half-width half-maximum (HWHM) of 0.14 ± 0.04°. For 
nuclear Bragg reflections from the crystal, the HWHM is 
resolution-limited at 0.155 ± 0.001°, indicating that the 
magnetic peak is resolution limited. This means that the 
coherence length of the magnetic order is greater than 
2400 A along the rock direction, which is at an angle of 
50° to the field. 

Using all of the blades and the results for different rock 
angles, a scattering map can be constructed around the 
magnetic peak. The blades cover a region of the scat- 
tering plane around the central spot of approximately 
±(5/5, 5/5, 5), where 5 = 0.05 rlu. Because of the verti- 
cal acceptance of the blades, a region out of the scatter- 
ing plane of approximately ± (5,-5, 0) is also covered. 
The scattering map for the (110) nuclear peak gives the 
resolution ellipsoid, which has its longest dimension ap- 
proximately parallel to the total momentum transfer Q. 
A fit to the magnetic peak shows that within errors its 
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Figure 3: (color online) Scan of the sample rotation angle uj 
through the magnetic peak at (0.56 -0.44 1.5), rotating the 
sample relative to the field and the incident neutron beam, 
with the magnetic field close to 12° to the basal plane. The 
background is comprised of data taken at 10.9 T (all blades), 
and those blades showing no peak at 10.8 T. The horizontal 
line is a straight line fit to the background data, and the other 
line is a Gaussian fit to the peak on top of the horizontal 
background. The temperature of the sample was 50 mK. A 
monitor count of 8000000 represents approximately two hours 
of counting time. 



shape is resolution limited. There is no evidence in the 
region studied for an additional FFLO-induced modula- 
tion parallel to the field direction [l7J , although it should 
be emphasized that our signal is weak, and satellites due 
to any modulation are expected to be weaker still. 

One can evaluate the magnetic moment per Ce ion by 
comparing the intensity of this peak with that of the very 
weak (110) nuclear peak, assuming that it does not suf- 
fer from extinction. If extinction is important, the value 
obtained will represent an upper limit on the magnetic 
moment. Due to the geometry constraints, only two nu- 
clear Bragg peaks were accessible, the (1 1 1) and the 
(1 1 0). On comparing the intensity ratios of these two 
peaks to calculations of the expected intensity (including 
absorption corrections) , there is clearly extinction of the 
stronger reflection, the (1 1 1). However, the extent of 
the extinction of (1 10) could not be confirmed due to 
the experimental constraints. The paramagnetic contri- 
bution to the peak intensity of these nuclear peaks was 
not taken into account in the calculations below, but es- 
timates indicated that it would have a negligible (< 1%) 
effect. 



For the model presented in Ref. [13J, which has the mod- 
ulated moments along the easy c axis, our upper limit on 
the Ce moment is 0.17 ± 0.02 ^s/Ce ion at 10.8 T, with 
the field 12° from the basal plane. It is possible that, 
when the field is not in the basal plane, the modulated 
magnetic moments rotate to move perpendicular to the 



Figure 4: (color online) The field dependence of the magnetic 
peak (0.56 -0.44 1.5) with a field angle of 12°. The back- 
ground is the averaged sum of the data obtained from the 
other analyzer blades during the measurement, and so rep- 
resent neighboring momentum transfer values. Corrections 
for blade efficiency have been made. The horizontal line is a 
straight line fit to the background points. The vertical lines 
mark the transition fields as observed from magnetostriction 
measurements in Ref. fig. 



applied field. We find that this is a small effect for a 12° 
tilt; the upper limit on the moment per Ce ion would be- 
come 0.16 ± 0.02 hb- These values are comparable with 
those found by Kenzelmann et al. [l3, 14 1, but because of 
the large absorption and uncertainty due to extinction, 
we hesitate to compare absolute magnitudes. 

The field dependence of the peak was also measured 
(Figured]), by sitting on the top of the peak and altering 
the field. Limited rocking scans were carried out at 10.9 
T and 9.5 T to check that the disappearance of the peak 
intensity was not due to its position moving with field. 
The field-dependence of the peak maps very well onto the 
phase diagram mapped out by magnetostriction (lp| . as 
do the results of Kenzelmann et al. [l3| • 

The other field angle we studied was 17°, over a field 
range from 8.5 T to 10.5 T, and at a temperature of 50 
mK. No magnetic peak was observed, probably because 
the magnetic moment has declined with field angle such 
that there is insufficient contrast. Using the approach 
described above, an upper limit of 0.025 (1b can be cal- 
culated for the modulated moment. It is unlikely that 
the loss of intensity is due to a shift of the peak posi- 
tion, because no change in the magnetic wavevector is 
observed between 0° and 12°. Due to the finite volume 
of reciprocal space surveyed, we would detect any small 
shift between 12° and 17°. 

In an antiferromagnetic structure, we expect that the 
modulated magnetic moments will lie along the easy axis. 
If a magnetic field is applied, the modulated moments will 
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tend to align perpendicular to the field. The magnetic 
structures reported in Refs. 13 and [13 clearly satisfy both 
of these conditions: in both cases the field is in the basal 
plane and the moments point along the c axis. However, 
when the magnetic field is rotated towards the c axis, as 
is done here, these conditions cannot be met simultane- 
ously, and we expect the magnetic order to become less 
stable. 



Our results, and the field dependent data in Ref. 113 . 
both match up with the phase diagram deduced from 
the magnetostriction measurements of Correa et al. [Io| , 
indicating that the magnetostriction signature correlates 
with the existence of the magnetically ordered Q-phasc. 
This signature weakens and has disappeared when the 
field is applied more than 22° out of the basal plane, 
suggesting that this is the limiting field angle for the oc- 
currence of the Q-phase. From our limited measurements 
at higher angles, we find no evidence for magnetic order- 
ing. These two experiments together indicate that the 
Q-phase ordering gets weaker and disappears as the field 
is rotated out of the basal plane. We therefore conclude 
that the Q-phase is not related to the anomalies seen in 
numerous bulk measurements when the field is applied 
parallel to the c axis; the phase associated with these 
anomalies must have another origin, and may be a true 
FFLO phase. 

The origin of the magnetic order in the Q-phase re- 
mains an open and very interesting question. In this 
phase, the ordering wavevector q mag = (0.44,0.44,0.5) 
is effectively unchanged on rotating the applied field di- 
rection - both in and out of the basal plane, and is close 
to that observed when superconductivity is suppressed 
by Cd-doping Jl8j . A spin resonance peak is seen at 
QS-R = (jii) [12|; and moves to lower energy as a field 
is applied [20( . These all point towards the ordering being 
controlled by a peak in susceptibility at a characteristic 
Fermi surface nesting vector [2l[, presumably between 



the dominant quasicylindrical sheets running parallel to 
c* 22|. Nonetheless, the appearance of the Q-P nase only 



within the superconducting state has to be explained. 
The repeat of the magnetic structure in the basal plane 
is a/( yf(2)q) ~ 7 A. The distance between planes of min- 
imum Ce moment in the Q-phase = o/(v2x 2(0.5 — q)) ~ 
27 A. Both of these lengths are considerably smaller than 
the superconducting coherence length, indicating that 
the dominant order parameter is prob ably not spatially- 
varying triplet superconductivity [U 
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Instead, the 

disappearance of the Q-phase in the normal state, and 
its suppression when the field is moved out of the basal 
plane suggests that the magnetic order just becomes sta- 
ble because of a small enhancement of the antiferromag- 



netic susceptibility in the mixed state [21|, that may be 



as the field direction is moved towards c. The confirma- 
tion of this model requires greater intensity to observe 
satellite magnetic peaks expected parallel to the applied 
field, or NMR measurements as a function of field angle. 
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brought about by strong Pauli paramagnetic depairing 
251 ] . or is perhaps driven by an FFLO-type pair density 
wave [l7l 12(1 ] , which would survive the magnetic ordering 



[1] 

[2] 
[3] 
W 
[5 
[6 
[' 
[8 
[9 
[10 

[11 
[12 

[1.3 
[14 

[15 
[16 
[17 

[18 
[19 
[20 
[21 
[22 
[23 

[24 

[25 

[26 



C. Petrovic et al, J. Phys.: Condens. Matter 13, L337 
(2001). 

R. J. Ormeno et al, Phys. Rev. Lett. 88, 047005 (2002). 
Y. Kasahara et al, Phys. Rev. B 72, 214515 (2005). 
Izawa et al, Phys. Rev. Lett. 87, 057002 (2001). 
Kumagai et al, Phys. Rev. Lett. 97, 227002 (2006). 
A. Radovan et al, Nature 425, 51 (2003) 
Bianchi et al, Phys. Rev. Lett. 91, 187004 (2003). 
Watanabe et al, Phys. Rev. B 70, 020506(R) (2004). 
Capan et al, Phys. Rev. B 70, 134513 (2004). 
F. Correa et al, Phys. Rev. Lett. 98, 087001 (2007). 
Fulde and R. A. Ferrell, Phys. Rev. 135, A550 (1964). 
I. Larkin and Yu. N. Ovchinnikov, Soviet Physics 
JETP 20, 762 (1965). 

M. Kenzelmann et al, Science 321, 1652 (2008). 

M. Kenzelmann et al, Phys. Rev. Lett. 104, 127001 

(2010). 

N. J. Curro et al, J. Low Temp. Phys. 158, 635 (2010). 
K. Lefmann et al, Physica B 385-386, 1083 (2006). 
Y. Yanase and M. Sigrist, J. Phys. Soc. Jpn. 78, 114715 

(2009) . 

M. Nicklas et al, Phys. Rev. B 76, 052401 (2007) 

C. Stock et al, Phys. Rev. Lett. 100, 087001 (2008). 
J. Panarin et al. J. Phys. Soc. Jpn. 78, 113706 (2009). 
K. Machida, private communication (2010). 

A. Koitzsch et al, Phys. Rev. B 79, 075104 (2009). 

A. Aperis et al, J. Phys. Condens. Matter 20, 434235 

(2008). 

D. F. Agterberg, M. Sigrist, and H. Tsunetsugu, Phys. 
Rev. Lett. 102, 207004 (2009). 

R. Ikeda, Y. Hatakeyama, and K. Aoyama, 
larXiv:1003.0309V 2 (2010). 

G. Koutroulakis et al, Phys. Rev. Lett. 104, 087001 

(2010) . 



